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Abstract 

We present a unified picture of mesons and baryons in the Dyson-Schwinger/Bethe-Salpeter approach, wherein the 
quark-gluon and quark-(anti)quark interaction follow from a systematic truncation of the QCD effective action and 
includes all its tensor structures. The masses of some of the ground state mesons and baryons are found to be in 
reasonable agreement with the expectations of a ‘quark-core calculation’, suggesting a partial insensitivity to the details 
of the quark-gluon interaction. However, discrepancies remain in the meson sector, and for excited baryons, that suggest 
higher order corrections are relevant and should be investigated following the methods outlined herein. 
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1. Introduction 

Hadrons provide a rich experimental environment for 
the study of the strong interaction, from details of the 
resonance spectrum to form factors and transition decays 
via electromagnetic probes. These reflect the underlying 
substructure of bound states by resolving, in a non-trivial 
way, the quarks and gluons of which they are composed. 
A theoretical understanding of hadrons in terms of these 
underlying degrees of freedom, interacting as dictated by 
quantum chromodynamics (QCD), is an on-going effort. 
Many approaches tackle it in different ways, simplifying 
certain aspects of the theory. Probing sensibly our the¬ 
oretical constructs with experimental input thus provides 
understanding of the theory itself. 

In continuum approaches to QCD, it is not possible in 
general to include all possible correlation functions in a 
calculation, as there are infinitely many of them. Al¬ 
though this can be viewed as a limitation of continuum 
approaches, only a finite number of these correlation func¬ 
tions have a significant role in the observable properties 
of hadrons. Therefore, by including a greater number 
of relevant correlation functions into the system, contin¬ 
uum methods provide an ideal framework to unravel the 
underlying mechanisms that generate observable effects 
from the elementary and non-observable degrees of free¬ 
dom of QCD. This is in contrast to lattice QCD calcula¬ 
tions, which can be viewed as theoretical experiments in 
the sense that, although they contain a priori all the dy¬ 
namics of QCD, it is challenging to single out individual 
contributions to a particular measurement. This makes 
these two approaches complementary. 
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Amongst the different continuum approaches, the com¬ 
bination of Dyson-Schwinger (DSE) and Bethe-Salpeter 
equations (BSE) has proven to be extremely useful in the 
calculation of hadronic properties from QCD [IH2- Typ- 
ically, solutions of DSEs constitute the building blocks 
(propagators and vertices) of bound-state calculations us¬ 
ing BSEs, which provide the bridge between QCD and 
observables. As described in more detail below, the inter¬ 
action terms that are kept in the DSE determine the in¬ 
teraction kernels among constituents in the BSEs, thereby 
defining a particular truncation of the DSE/BSE system. 
One works towards a model independent truncation by 
including a larger set of interaction terms; although this 
programme is obviously not achievable in its totality, it is 
expected that there will be some degree of convergence at 
the level of this vertex expansion. 

Here we focus upon the inclusion of the quark-gluon 
interaction, the reliable construction of which is a chal¬ 
lenging task. However, one is guided by various symme¬ 
tries - notably that of chiral symmetry - that provide for 
stringent constraints. To implement these symmetries at 
the level of the quark and gluon interaction, simplifica¬ 
tions are clearly necessary which typically fall into three 
categories: (i) the quark-gluon vertex is truncated to its 
tree-level component times a momentum-dependent effec¬ 
tive coupling with the quark DSE and hadron BSEs solved 
self-consistently mu-, (ii) a more sophisticated model for 
the quark-gluon vertex is used, with the contribution from 
its different tensor structures modelled, hence abandon¬ 
ing self-consistency but gaining instead flexible insight into 
the relative importance of each of these structures [Mai; 
(iii) some non-perturbative effects of the quark and gluon 
interaction are taken into account by solving the quark- 
gluon vertex DSE self-consistently, but potentially intro¬ 
ducing some truncation artifacts [TOUT^ . We follow here 
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the latter approach, since as we demonstrate in this letter 
it enables the controlled inclusion of interaction mecha¬ 
nisms based on a loop expansion of the effective action. 

The effective action is a generating functional for proper 
vertex functions, and may be considered as a means to de¬ 
fine the quantum field theory given an action, since all 
necessary Green’s functions can be derived from it. Re¬ 
lated are the n-particle irreducible (nPI) effective actions 
that form a family of different, but equivalent, represen¬ 
tations of the same generating functional (see, e.g. USD- 
They are defined as functionals of all m < n Green’s func¬ 
tions of the theory (fields, propagators, vertices, etc.). Al¬ 
though its exact form is not known in general, its loop 
expansion in h is well-defined m and can in practice be 
performed. Moreover, each term in the expansion already 
captures both perturbative and non-perturbative physics. 

One reason that makes nPI techniques a powerful tool 
is that they provide a natural link between bound-state 
equations described by BSEs and the propagators and ver¬ 
tices provided by DSEs pAhTS] . A truncation of the loop 
expansion at a certain order translates into a unique pre¬ 
scription for the truncation of the DSEs and the BSEs 
that maintains symmetries. Eor the study of two- and 
three-body states, it suffices to use either the 2PI effective 
action, which is defined in terms of fully-dressed propa¬ 
gators but bare vertices m, or the 3PI effective action, 
which is defined also in terms of the fully-dressed vertices. 
In this work we restrict ourselves to the 2PI case and de¬ 
fer the use of the 3PI effective action to a future and more 
comprehensive study. 

In connection with the three categories outlined above, 
it is worth mentioning here that the somewhat hybrid pos¬ 
sibility of supplementing some of QGDs degrees of freedom 
in favor of effective ones, such as pions, has also been ex¬ 
plored [2Qll25] . These can be viewed as approximate rep¬ 
resentations of the four-quark vertex in the 4PI formalism 
that introduces at the first step decay channels and a mix¬ 
ing with tetraquark states. 

In the present work, we incorporate the results of a 
recent study of the quark-gluon vertex from a truncated 
DSE [12] in the calculation of meson and baryon masses. 
That truncation can be interpreted in the context of the 
2PI effective action at 3-loop. Although on the technical 
side this is no novelty for meson calculations miEn], it 
is the first time that corrections incorporating the gluon 
self-interaction have have been included in the covariant 
three-body baryon calculation. While there exist other re¬ 
cent investigations of the quark-gluon vertex [?7ll29] , these 
have not yet been confronted with the challenge of repro¬ 
ducing hadron phenomena for reasons we discuss below. 

Einally, we wish to stress that this work represents only 
the first step in an on-going effort to incorporate realistic 
QGD’s Green’s functions into the self-consistent study of 
bound states. Not surprisingly, the low-order of the trun¬ 
cation used performs only as well as simple phenomeno¬ 
logical models such as rainbow-ladder. However, but most 
importantly, it serves as a proof of principle that such an 



Figure 1: The 2-particle irreducible term r 2 [T,G] in the definition 
of the effective action, up to three loops. 
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Figure 2: The Dyson-Schwinger equation for the quark propagator. 


endeavor is feasible, as further increasing the order of the 
truncation does not increase the technical complexity dra¬ 
matically. 


2. Framework 

The starting point for the study of hadronic observables 
in the present framework is thus the effective action 

r[«',G] =S[^']+i'IVlogG-i'IVGo^G + r2[4',G] , (1) 

where S is the classical action, and T and G collectively 
represent the fields and full propagators of QGD, respec¬ 
tively. The term r 2 [T, G] contains two-particle irreducible 
diagrams only and Go denotes the classical propagators. 
To proceed, we perform a loop expansion of r 2 [T,G] to 
three-loop order, as shown in Eig. Moreover, we keep 
only the non-Abelian term that connects the gauge to the 
matter sector and neglect the Abelian correction (third di¬ 
agram), as it is expected to be subleading in the large-A^c 
limit; whether this is indeed the case for the description 
of hadron phenomena must certainly be tested and will be 
the subject of future work. 

The next basic element in meson and baryon calcula¬ 
tions is the fully-dressed quark propagator. It is given by 
the quark DSE, see Eig. 

S-\p-,^i) = Z2So\p) + ^{p;fi) , (2) 

with quark self-energy 

Sfe p) = g^Z.fCp [ rS{q)r^{q,p)D^,(k) . (3) 

Jk 

Here g = /c +p, the integral measure is = f d^k}(27r)^ 
and Zif, Z 2 are renormalization constants for the quark- 
gluon vertex and quark propagator respectively. It is 
clearly dependent upon both the gluon propagator 
and the quark-gluon vertex r^(g,p). The (Landau gauge) 
propagators are 

S~'^{p) = zj^{p^) {if + Mip"^)) , (4) 

= Ti;^^z{e)ie , (5) 
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with quark wavefunction Z^ dynamical mass M{p^) 

and gluon dressing Z{k‘^). The transverse projector is 
"^(k) ~ ~ From the 2PI effective action, the 

quark DSE is determined via a functional derivative with 
respect to the quark propagator 


E = 



(6) 


The expansion of the effective action in Fig. [T] thus defines 
a truncation of the quark DSE, which is equivalent to the 
truncation of the quark-gluon vertex DSE shown in Eig.[^ 
Specifically, the truncated vertex DSE can be given as a 
summation of vertex corrections 


r^il,k) = Z^f^^+A^^ + ... , (7) 

with the non-Abelian term (see Eig. |^, and the el¬ 
lipsis denotes contributions not considered here. The non- 
Abelian correction has the form 

N f- - 

Ana(^. ^) = T y (/2, -q2) (8) 

xrfg^'^{qi,q2,P3)D"^' {qi)D^^' {q2) , 


where the internal vertices T^ retain only the tree-level 
tensor structure and are supplemented by an enhancement 
factor in the infrared. Such an enhancement is needed in 
order to account for the effect on the quark-gluon vertex of 
higher loop terms in the 2PI effective action. An alterna¬ 
tive to this is using a 3PI or higher effective action to define 
the truncation, which means that both vertices and prop¬ 
agators are fully dressed and independent objects [ISlIIHl. 
Eurther details on ^ and its solution can be found in [12] . 

A second functional derivative of the effective action 
with respect to the quark propagator defines the quark- 
antiquark BSE interaction kernel [HHini na na. 

Equivalently, 




J ik-,lj 




[E {x',y')] 


(9) 


followed by a Eourier transform to momentum space. This 
is the kernel that appears in the BSE description of a me¬ 
son as a quark-antiquark bound state 


[rM(p,P)]« = f , (10) 

J k 


where Vm{PiP) is the Bethe-Salpeter amplitude and 
XM{k,P) = S{k-^)TM{k^ P)S{k-) its wavefunction. The 
quantum numbers of the state are defined by its covari¬ 
ant decomposition [3Qf[32] . This gives access to both the 
mass of the bound-state as well as details of its internal 
structure. 

The crucial motivation for such a definition of the 
Bethe-Salpeter kernel is that chiral symmetry, as ex¬ 
pressed via the axial-vector Ward-Takahashi identity 



Figure 3: The truncated DSE for the quark-gluon vertex, showing 
the non-Abelian correction. Boxes represent an RG improvement of 
the bare vertex; see 



Figure 4: The quark-antiquark kernel beyond rainbow-ladder. 


(axWTI) 

[5^(p+) 75 + 75S(p-)] = (11) 

J k 

is correctly implemented in the calculation of meson prop¬ 
erties [33]. This guarantees, in particular, the identifica¬ 
tion of the pion as a Goldstone boson in the chiral limit. 

With the truncation of the vertex DSE given in (|^, the 
cutting process gives rise to a quark-(anti)quark kernel, 
given in Eig. whose topology is analogous to a gluon 
ladder 


( 12 ) 

This, taken with the truncated quark-gluon vertex above, 
satisfies chiral symmetry by construction. Eurther correc¬ 
tions to the quark-gluon vertex, in particular those featur¬ 
ing additional quark propagators, yield different topologies 
and/or are higher loop order; they are beyond the scope 
of the present work but can in principle be included. 

It is important to stress that when the 2PI effective ac¬ 
tion is used to define the truncation at the level of ver¬ 
tex functions, the expansion must be chosen such that 
the functional derivative can be formally performed to 
define the chiral-symmetry preserving BSE kernels. This 
implies that it must be possible to resolve diagrammati- 
cally the quark lines in the vertex corrections 0 . This is 
the reason why quark-gluon vertices defined as ansatze are 
difficult to accommodate (and, hence, test) consistently in 
this framework. 




3 











It is illustrative to consider here the case of the two-loop 
expansion of the 2PI effective action, which leads to the 
well-known rainbow-ladder truncation. There, suppressing 
constants and color factors for simplicity, we have 

[I](p)] C. [ [j'^S{q)r]..D^.{k) , (13) 

J k 

from which a functional derivative provides the kernel 

( 14 ) 

In practical calculations the gluon propagator is modelled 
by a phenomenological function which includes the afore¬ 
mentioned infrared enhancement of the quark-gluon vertex 
from (omitted) higher order corrections. 

Since we wish to also consider bound-states of three- 
quarks, it proves useful as an intermediate step to formu¬ 
late the diquark bound-state 


[rc(p,P)] = f [xD(fe.P)]„ , (15) 

J k 

with diquark amplitude T ^ (p, P) and wavefunction 
Xoik^P) = S{k-^)r]j{k, P)S'^{—k-). The superscript T 
denotes transposition; the covariant decomposition of a 
J~^ diquark Tq are that of a meson TmC together 
with a charge conjugation matrix. The corresponding di¬ 
quark kernel is then 


that is, it amounts to a transposition of the lower spin- 
line. For general ladder-like kernels, such as that in (14) 
the color factors are — l)/2A^c and —{Nc + l)/2A^c foi* 
a meson and diquark, respectively. 

The analogue of the Bethe-Salpeter equation for 
baryons can now be formulated. It contains the permuted 
sum of the two-body quark-quark kernel and an ir¬ 
reducible three-body kernel [151 IS] 


3 

^ ^-1 Gr(3)^ ^ (17) 

a=l 

see Fig-i We use here a compact notation, omitting in¬ 
dices, where implied discrete and continuous variables are 
summed or integrated over, respectively. Here Gq ^ repre¬ 
sents the product of three fully-dressed quark propagators 
S. The subscript a labels the quark spectator to the two- 
body interaction. At the truncation level of the 2PI effec¬ 
tive action we will consider in this work, the irreducible 
three-body force is identically zero [18]. The two-body 
kernel follows, as before, from the quark and quark-gluon 
vertex DSE. 


3. Results 

A consistency check that chiral symmetry is correctly 
implemented is provided by calculating the pion mass in 




Figure 5: The covariant Faddeev equation for the baryon, see |T7|. 


the chiral limit. Finding a massless pion, as we indeed do, 
serves as a verification that the numerics are under control. 

In [T2| the internal vertices were adjusted to both resem¬ 
ble the tree-level component of the calculated vertex, and 
to yield a running quark mass in agreement with lattice 
calculations. Once the vertex truncation is defined, the 
only free parameter left is the current-quark mass. The 
light quark mass is then fixed to rriu/d = 3.7 MeV at the 
renormalization point 19 GeV so that the physical pion 
mass is obtained. 

In Table [^ we show the ground-state meson masses be¬ 
low 1.4 GeV for two quark flavors for rainbow-ladder (us¬ 
ing the Maris-Tandy interaction 0) and for the beyond 
rainbow-ladder truncation presented here. The case of 
the scalar and vector states is remarkable in that their 
masses agree well with the experimental value of the p- 
and ( 7 - mesons; this is, however, an unexpected feature. 
On the one hand, it is well known that unquenching ef¬ 
fects (absent in this calculation) such as pion-cloud cor¬ 
rections ED mi Eg and effects associated with decay pro¬ 
cesses ED EH provide extra attraction for vector and/or 
scalar channels. On the other hand, it is not yet settled 
whether the physical cr-meson originates (at least in part) 
form a qq state. Moreover, even if that was the case, 
the measured mass would be the result of the dynami¬ 
cal coupling of this hare state with several other channels 
[391140]. In this respect, it is interesting to note that sev¬ 
eral DSE/BSE calculations [6l|9l|32], including the present 
one, generate a relatively light scalar qq state and hence 
suggest that the picture of a purely dynamical origin of 
the cr-resonance is not realised; this can likely be traced to 
the back-coupling of the quark-gluon vertex to the vertex 


Table 1: Meson masses and pion decay constant in GeV as calculated 
in rainbow-ladder [5] and beyond rainbow-ladder. Results affixed 
with are fitted values. _ 



RL 

BRL 

PDG ES 

0 (tt) 

0++ (a) 
1-- (P) 

1++ (ai) 

1+- (6i) 

0.14t 

0.64 

0.74 

0.97 

0.85 

0.14t 

0.52 

0.77 

0.96 

1.1 

0.14 

0.48(8) 

0.78 

1.23(4) 

1.23 

U 

0.092T 

0.103 

0.092 
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Table 2: Baryon masses in GeV as calculated in rainbow-ladder El 
mss] and beyond rainbow-ladder. 



RL 

BRL 

PDG [23 

1 / 2 + (N) 

0.94 

1.05 

0.94 

1 / 2 - (N*) 

1.16 

1.08 

1.54(1) 

3 / 2 + (A) 

1.22 

1.24 

1.23 


itself and will be remedied at higher order in loops or in 
nPL 

The previous remarks can be summarized in that one 
would expect our quark-core calculation to lead, in general, 
to heavier meson masses. The fact that this is not the 
case is an indication that other mechanisms, absent in the 
present calculation, are of relevance even in the ground- 
state meson sector. Additional support for this is found 
in the ai/bi — p and ai — bi splitting; whereas the former 
is improved with respect to the RL value, the latter is 
exceedingly large showing an imbalance in the different 
components of the present quark-gluon vertex. 

Let us now turn our attention to baryons. In Table 
we show the calculated nucleon and delta-baryon masses. 
We also calculated the mass of the nucleon parity partner 
1 / 2 “, as it is the first signature of the failure of RL in 
the baryon sector. Clearly, the masses for positive-parity 
states are slightly overestimated. This suggests that the 
simple quark-core picture supplemented with pion-cloud 
effects [23l [24l |3^ |44] could be realized for ground-state 
baryons. The possibility that ground-state baryons are, 
in contrast to mesons, less sensitive to the details of the 
quark-gluon interaction is an interesting one. It explains 
why the simple rainbow-ladder truncation has been thus 
far so successful in describing baryon properties [Ml SU¬ 
SS]. 

The situation for excited baryons does not appear to 
be so simple. In RL the negative-parity nucleon comes 
out extremely light [42] unless one inflates the interaction 
strength |45l SH], thus spoiling the agreement for other 
observables. The present calculation, using a genuine (al¬ 
beit truncated) solution of a QCD-vertex DSE, does not 
improve on the situation of the nucleon’s negative-parity 
partner. This result by itself already makes a case for 
the systematic introduction of interaction mechanisms, as 
it is not known a priori whether missing Abelian and/or 
non-Abelian mechanisms are leading in determining the 
spectrum of (bare) baryon resonances. 

Another interesting aspect of using as interaction a solu¬ 
tion of the quark-gluon vertex DSE is that it naturally fea¬ 
tures a quark-mass dependence of the interaction strength. 
One class of observables in which such a mass dependence 
is manifest are the baryon sigma-terms. Calculated here 
using the Eeynman-Hellmann theorem, they measure the 
dependence of the baryon mass on the quark mass. As 
shown in Tablethese are exceedingly small in RL, owing 
to the quark-mass independence of the effective interac¬ 
tion. This changes dramatically in the present calculation. 


Table 3: Sigma terms in MeV as calculated in rainbow-ladder and 
beyond rainbow-ladder. 



RL 

BRL 

Other [ 47 II 55 ] 

1 / 2 + 

(AT) 

30 

60 

25-60 

3 / 2 + 

(A) 

24 

63 

32-79 


Eor the nucleon we now obtain a value in reasonable agree¬ 
ment with the upper limit of the consensual range HMI]. 
Eor the delta, although there is no well established value to 
serve as a reference, we find a similarly large result for the 
sigma term. However, this could be related to the absence 
of decay mechanisms that induce non-analytic behaviours 
in the baryon-mass curve near the physical point. 

4. Summary 

In summary, we presented the first calculation of both 
mesons and baryons in a mutually consistent truncation 
of the quark-gluon vertex beyond rainbow-ladder. Being 
fixed in accordance to the 2PI effective action at three- 
loops, the freedom to adjust model parameters is limited 
and hence the bulk features of the results for the meson 
and baryon masses cannot presently be altered. Exten¬ 
sions must therefore be made that either take into account 
higher loop corrections to the effective action, or as we 
envisage, to incorporate vertices dynamically into the ex¬ 
pansion by considering 3PI actions and beyond m- 

Regardless, the means to connect quarks and gluons pre¬ 
sented here via symmetry preserving and calculated inter¬ 
action vertices is both an important technical and con¬ 
ceptual achievement that, being in itself independent of 
the particular ansatz, can be systematically applied to 
e.g. crossed ladder and other higher order corrections to 
the kernel. Being intrinsically dependent upon the quark 
flavor and hence the quark mass, it will form the basis 
of future investigations of the baryon octet-decuplet [43] 
and meson nonet for the strange quark [32] |56], as 
well as heavy-quark studies of charmonium and bottomo- 
nium [ 57 II 59 ] . Eurthermore, the impact of corrections 
beyond rainbow-ladder on the internal structure of the 
hadrons can be tested through the calculation of electro¬ 
magnetic form-factors [6Qll62] . 
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